Scarring in a driven system with wave chaos 
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We consider acoustic wave propagation in a model of a deep ocean acoustic waveguide with a 
periodic range-dependence. Formally, the wave field is described by the Schrodinger equation with 
a time-dependent Hamiltonian. Using methods borrowed from the quantum chaos theory it is shown 
that in the driven system under consideration there exists a "scarring" effect similar to that observed 
in autonomous quantum systems. 



It is well-known in the theory of quantum chaos that 
classical periodic orbits, to a significant extent, determine 
both the distribution of energy levels and structures of 
eigenfunctions [llliy . In particular, for some eigenfunc- 
tions, the amplitudes are high in the vicinity of unstable 
periodic orbits and low elsewhere. This phenomenon, 
first discovered in a quantum billiard, is called the scar- 
ring Q. It has been established that the scarring is a 
common property for many autonomous chaotic systems 
0, 0, 0. The description of scarring is related to the 
problem of construction of eigenfunctions on the basis of 
purely semiclassical calculations 0, H S • 



obeys the parabolic equation [2y, [2 
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It is also well-known that so different physical problems 
as quantum particle dynamics and wave propagation in 
nonuniform media in the narrow beam approximation are 
described by the same parabolic (or Shrodinger) equa- 
tion. One of the most intriguing object is the underwater 
acoustic (UWA) and the long range sound propagation in 
the ocean 0, [3 • Application of methods of quan- 
tum chaos to the wave chaos proves to be efficient in 
many publications E E E [la 20]. The goal of this 
paper is to extend this analogy between quantum and 
wave chaos and to present simulations that demonstrate 
the phenomenon of scarring in a nonautonomous Hamil- 
tonian system with 1.5 degrees of freedom related to the 
UWA. We consider a simplified case of wave propagation 
in a two dimensional range-dependent model of hydroa- 
coustic waveguide that has been used to study wave and 
ray chaos at lon g ra nge sound propagation in the ocean 
El 111 111 in |2|. Formally, the simplified model of 
ray propagation is equivalent to a nonlinear oscillator 
perturbed by a periodic force, while the corresponding 
parabolic equation is equivalent to the Shrodinger equa- 
tion for the oscillator. 



Consider a two-dimensional UWA waveguide with the 
sound speed c as a function of depth, and range, r. 
In the parabolic equation approximation (valid under as- 
sumptions that waves propagate at small grazing angles) 
the monochromatic wave field at a carrier frequency / 



where U{r,z) = (l — c^/c^(r, z)) /2. Here c and k = 
27r//c are a reference sound speed and a reference 
wave number, respectively. Equation Q formally coin- 
cides with the time-dependent Schrodinger equation. In 
this analogy r plays a role of time, H is the Hamilto- 
nian operator, U{r^z) is the potential, and k~^ asso- 
ciates with the Planck constant. The ray paths are de- 
termined by the standard Hamilton equations dp/dt = 
—dH/dz and dz/dr = dH/dp with the Hamiltonian 
H = p'^/2 + U{r^z)^ where the momentum p is related 
to the ray grazing angle hy p = tan 0. 

Following Ref. 21] we take the sound speed field 
c(r, = co{z) +(5c(r, z), where co{z) is the so-called 
Munk profile widely used to model sound propagation 
in the deep sea 27], 

co{z) = c (l + £( (e-^(^-^-)/^ + 2{z - Za)/B - 1 

"(2) 

with c = 1.5 km/s, Za = —1 km, e = 0.0057, and B = 1 
km, and 

Sc{r,z) = -2-fcz/Be-^^^^ sin(27rr/L). (3) 

with 7 = 0.01 and L = 10 km is a perturbation that 
generates the ray chaos. The corresponding potential 
U (r, z) is a periodic function of range. Its shapes at three 
different ranges are shown in Fig. 1. 

Figure 2 presents the phase portrait of the ray sys- 
tem plotted using the Poincare map ^1 0, 0, HI HI] 
(pn+i,^n+i) = T{pn^Zn), whcrc T dcuotcs transforma- 
tion of the momentum and coordinate of a ray trajectory 
taken at range nL to range (n + 1)L for n = 0, 1, . . .. We 
see a typical picture of stable islands formed by regular 
curves submerged in a chaotic sea filled with randomly 
scattered points depicting chaotic rays. 

To find out how the structure in Fig. 2 reveals itself 
in a wave dynamics described by Eq. (QJ we shall use 
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FIG. 1: Unperturbed potential U{0,z) (solid line) and its 
maximal deviations (dashed line) due to the perturbation. 
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FIG. 2: Phase portrait of the ray system. 

the Floquet theorem 0, and follow Ref. [iOI where a 
similar approach was applied. The solution of Eq. (QJ 
can be presented as a sum of Floquet modes 



,(r,z) = e^"'»'-/^*™(r,^), 



(4) 



where ^rn(^ L^z) = ^^(r, z) are range-periodic, am 
are real constants, and m = 0, 1, . . .. In quantum me- 
chanics analogous terms describe the so-called Floquet 
states, and constants a^i are proportional to the quasi- 
energies. ^rn(^7^) are eigenfunctions of the shift opera- 
tor F: F^rnir, z) = ^^(r + L, 2:) = e'^'^^'^ rnir, z). Since 
F is Hermitian, ^^(r, z) form a complete set 2, 3j. 

Let us consider functions il)m{z) = ^m(0, 2;). They can 
be presented as a decomposition tpm{z) = Cmq^qi^)^ 
where (Pq{z) {q = 0,1,...) are solutions to Sturm- 
Liouville eigenvalue problem Ho^Pm = Fm^m{z) in the 
unperturbed waveguide with the Hamiltonian Hq corre- 
sponding to c(r, z) = cq{z). The eigenfunctions (pm{z) 
and eigenvalues Em are readily found using a standard 
mode code [sj. It is easy to see that the phase fac- 
tors e*^"^ and columns of matrix ||C^g|| are eigenvalues 
and eigenvectors, respectively, of a unitary matrix with 



-I 



dz (pmiz)Fipq{z). 



(5) 



Note that F(pq{z) is a solution to the parabolic equa- 
tion at range r = L found for an initial condition 
u{0^z) = (fin (z). We have solved Eq. (QJ using the code 
MMPE Numerical results presented below have 

been obtained for a carrier frequency / = 200 Hz. 

Associating a parameter fi with each function V^rn(^), 
where /i = q \ Cmq\^ is sum of intensity weighted num- 
bers of eigenfunctions (pq^ we sort the Floquet modes in 
the order of increasing values of ja. Therefore the eigen- 
functions tl^m{z) with small m describe waves localized 
at depths z in the vicinity of the minimum of /7(0, z). 
The greater is the mode number m the steeper are the 
grazing angles of waves forming the Floquet mode. 

In order to link the structure of the Floquet mode with 
the classical phase space we shall consider the Husimi 
distribution function 0, IT^ 



Wm{p,z) 



X exp 



ikp{z' — z) 



dz'lpm{z') 



{z' - zf 
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(6) 



representing a projection of the wave field onto a Gaus- 
sian state with minimum uncertainty (coherent state) 
centered at a point (p, 2;). Since the ray grazing angle 
d with respect to the r-axis is related to the momentum 
by ^ = arctanp the distribution Wm{p^ z) may be inter- 
preted as a local spectrum of the wave field smoothed 
over spatial and angular scales. 

Figure 3 shows Husimi functions evaluated for six 
Floquet modes. These examples have been selected to 
demonstrate a variety of structures of Floquet modes in 
our environmental model. At each panel we have out- 
lined borders of stable islands discernible in the central 
part of the phase portrait. Although the range variations 
of the potential U are most pronounced at its minimum, 
the first Floquet eigenfunctions V^rn(^) are most closed to 
their unperturbed counterparts, i.e. to functions (fm{z) 
with the same numbers m. In the WKB approximation 
the Husimi function for any (pm{z) is a fuzzy version of a 
closed curve representing a period of an unperturbed ray 
whose Hamiltonian is equal to Em- In Fig. 3a we see an 
example of such a fuzzy curve. Here a Husimi function 
for V^io(^) is plotted. It almost coincides with the un- 
perturbed Husimi function for (pio{z) (not shown). The 
Husimi functions for the first 50 eigenfunctions V^rn(^) are 
localized within the central island. The fuzzy curves cor- 
responding to m approaching 50 repeat the shape of the 
border of the central island. This is clearly seen in Fig. 
3b where the Husimi function constructed for V^4o(^) is 
shown. 
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space visited by predominantly chaotic or regular rays, 
respectively. Moreover, Fig. 3c presents a Husimi dis- 
tribution localized outside the five sub-islands but in 
the vicinity of their borders. We conjecture that this 
is a finite wavelength analog of the stickiness to the 
borders of islands-surround-islands, i.e. the presence 
of such parts of the chaotic ray trajectory that cor- 
respond to the almost regular rotation of rays around 
the islands 34]. This occurs when, after wandering in 
the phase space, the trajectory approaches sub-islands 
and "sticks" to their borders for a fairly long time. 
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FIG. 3: Husimi distributions for the Floquet eigenfunctions 

'ipioiz) (a), ip4o{z) (b), ip66{z) (c), ip77{z) (d), ip88{z) (e) , and 
'ipiosiz) (f). Parameter Az = 0.01 km. 
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FIG. 4: Unstable (upper panel) and stable (lower panel) pe- 
riodic rays with the period 50 km. 



Figures 3b-f demonstrate that the Husimi functions 
may have quite different structures. But there exists 
a remarkable common feature typical for all the first 
140 Floquet modes considered in our numerical simula- 
tions: each distribution is localized within (i) the cen- 
tral stable island (Figs. 3a and b), (ii) within five sub- 
islands surrounding the central one (Fig. 3e), or (iii) 
within the chaotic sea (Figs. 3c, d, and f). Thus, the 
structure of Husimi functions correlates with the shape 
of stable islands. Roughly speaking, one may say that 
there are "chaotic" and "regular" Floquet modes whose 
Husimi functions are localized within areas of the phase 
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FIG. 5: Function |^77(r, z)| at a period of 
its range-dependence. Five segments of the 
unstable ray from the upper panel of Fig. 4 
are superimposed on the contour plot at the 
bottom. 

Figures 3d and 3e need more discussions. Each of them 
has five clearly seen maxima. In Fig. 3d the maxima 
are located within the chaotic sea, near the destroyed 
saddles, while in Fig. 3e they are located near the el- 
liptic points of the sub-islands. We argue that these 
maxima are associated with trajectories of two periodic 
rays graphed in Fig. 4. Both rays have periods equal 
to 5L = 50 km, while there are no other rays with the 
period of 50 km. The ray shown in the upper panel of 
Fig. 4 belongs to the chaotic sea. The difference \Sz\ 
between its coordinate and that of a ray with infinites- 
imally close starting parameters, on the average, grows 
as exp(Ar) with the Lyapunov exponent A = 1/60 km~^. 
In contrast, the ray in the lower panel has a zero Lya- 
punov exponent, i.e. it is regular. Both rays have been 
divided into five 10 km segments. The segments connect 
the consecutive waveguide cross-sections at range points 
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r = nL chosen for the stroboscopic observation of both 
the ray structure and the Floquet modes (at these range 
they are presented by ipmi^))- Smah black circles within 
red and yellow spots (corresponding to local maxima of 
Husimi function) in Fig. 3d indicate positions (in the 
phase plane) of starting points of five segments shown 
in the upper panel of Fig. 4. In agreement with our 
expectation these points lie at centers of five maxima of 
the Husimi function. The number of each segment is 
assigned to a corresponding circle and is indicated next 
to it in the figure. Similarly, the five maxima in Fig. 
3e are associated with starting points of segments of the 
stable ray shown in the lower panel of Fig. 4. It is nat- 
ural that the black circles depicting positions of these 
points are located in the centers of five stable sub-islands. 
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FIG. 6: The same as in Fig. 5 but for the 
Floquet mode ^88- Five segments superim- 
posed on the contour plot correspond to the 
stable ray from the lower panel of Fig. 4. 

The above results suggest that the two periodic rays 
shown in Fig. 4 determine essentially the structure of 
the Floquet modes with m =77 and 88. This conjec- 
ture is confirmed in Figs. 5 and 6 where we present 
distributions of intensities of these modes within a range 
interval (0,L), i.e. over a spatial period of all the Flo- 
quet modes. Both functions |^77(r, z)| and | ^88(^,^)1 , 
as well as | ^66(^,^)1 shown in Fig. 7 are scaled to have 
absolute maxima equal to one. In Fig. 5 we see two rep- 
resentations of function |^77(r, z)|. In the lower panel 
five segments of the unstable ray from the upper panel 
of Fig. 4 are superimposed on the contour plot repre- 
senting isolines of the intensity distribution. It is clearly 
seen that the intensity is high at the segments of the un- 



stable periodic ray and low elsewhere. So, the periodic 
unstable ray path scars the Floquet eigenfunction. This 
is an analog to the scarring phenomenon observed in dif- 
ferent autonomous quantum systems "i", (see also in 
(sj). In Fig. 6 it is seen that the 88-th Floquet mode 
localized within the five sub-islands also has relatively 
large amplitude in the vicinity of segments of the stable 
ray shown in the lower panel of Fig. 4. Moreover, prac- 
tically all functions |^m(^, z)\ with m > 50 display some 
pattern which may be associated with periodic rays sim- 
ilar to those shown in Figs. 5 and 6. We illustrate this 
statement in Fig. 7 for the Floquet mode with m = 66. 

Figures 5 and 6 suggest that the two Floquet modes 
may be considered as a superposition of five 10 km 
pieces of a wave beams propagating along the periodic 
unstable and stable rays, respectively. We conjecture 
that scars observed in intensity distributions of other 
modes, like that shown in Fig. 7, also may be as- 
sociated with segments of some periodic rays. How- 
ever, this issue needs a more detailed investigation. 
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FIG. 7: The same as in Fig. 5 but for the 
Floquet mode ^66 and without ray segments 
superimposed on the contour plot. 

The obtained results demonstrate that the scarring 
effect may be observed not only in autonomous quan- 
tum systems but also in driven ones described by the 
same equations as the waveguides with periodic range- 
dependencies. We presume that periodic rays determine 
the structures of the Floquet modes in approximately the 
same way as periodic orbits determine structures of the 
wave functions. The obtained results demonstrate that 
unstable chaotic rays can be considered as a reference di- 
rection for the corresponding wave packet propagation. 
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This effect may be used for the transfer of information 
and communication. 
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